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FIRST SEASON QUIET OBSERVATIONS: MEASUREMENTS OF CMB POLARIZATION POWER SPECTRA 

AT 43 GHZ IN THE MULTIPOLE RANGE 25 < i < 475 
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ABSTRACT 

The Q/U Imaging ExperimenT (QUIET) employs coherent receivers at 43 GHz and 95 GHz, oper- 
ating on the Chajnantor plateau in the Atacama Desert in Chile, to measure the anisotropy in the 
polarization of the CMB. QUIET primarily targets the B modes from primordial gravitational waves. 
The combination of these frequencies gives sensitivity to foreground contributions from diffuse Galac- 
tic synchrotron radiation. Between 2008 October and 2010 December, over 10,000 hours of data were 
collected, first with the 19-element 43-GHz array (3458 hours) and then with the 90-element 95-GHz 
array. Each array observes the same four fields, selected for low foregrounds, together covering pa 1000 
square degrees. This paper reports initial results from the 43-GHz receiver which has an array sensi- 
tivity to CMB fluctuations of 69 /iK-^/s. The data were extensively studied with a large suite of null 
tests before the power spectra, determined with two independent pipelines, were examined. Analysis 
choices, including data selection, were modified until the null tests passed. Cross correlating maps 
with different telescope pointings is used to eliminate a bias. This paper reports the EE, BB, and 
EB power spectra in the multipole range t = 25-475. With the exception of the lowest multipole bin 
for one of the fields, where a polarized foreground, consistent with Galactic synchrotron radiation, is 
detected with 3-<r significance, the E-mode spectrum is consistent with the ACDM model, confirming 
the only previous detection of the first acoustic peak. The B-mode spectrum is consistent with zero, 
leading to a measurement of the tensor-to-scalar ratio of r = 0.35+q 87- The combination of a new 
time-stream "double-demodulation" technique, side-fed Dragonian optics, natural sky rotation, and 
frequent boresight rotation leads to the lowest level of systematic contamination in the B-mode power 
so far reported, below the level of r = 0.1. 

Subject headings: cosmic background radiation — Cosmology: observations — Gravitational waves — 
inflation — Polarization 
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1. INTRODUCTION 

The inflationary paradigm resolves several outstand- 
ing issues in cosmology, including the flatness, horizon, 
and monopole problems, and it provides a compelling 
explanation for the origin of structure in the Universe 
(e.g. Liddle & Lyth 2000, and references therein). So far 
all cosmological data, including measurements of Cos- 
mic Microwave Background (CMB) anisotropies, sup- 
port this paradigm; still the underlying fundamental 
physics responsible for inflation is unknown. Inflation 
produces a stochastic background of gravity waves that 
induce odd-parity tensor "B modes" at large angular 
scales in the CMB polarization. If these primordial B 
modes, parametrized by the tensor-to-scalar ratio r, are 
detected, one can learn about the energy scale of infla- 
tion. In many attractive slow-roll models, this scale is 
given approximately by r 1 / 4 x 10 16 GeV. For large-field 
models, the energy scale is near the Grand Unification 
Scale in particle physics, so that r > 0.01. A new gen- 
eration of experiments aims for good sensitivity in this 
range of r. Establishing the existence of primordial B 
modes would both verify an important prediction of in- 
flation and provide access to physics at an incredibly high 
energy scale. 

The most stringent limit to date is r < 0.20 at the 95% 
confidence level (Komatsu et al. 2010) set by a combi- 
nation of CMB tcmpcrature-anisotropy measurements, 
baryon acoustic oscillations, and supernova observations, 
but cosmic variance prohibits improvements using only 
these measurements. 

E-mode polarization has now been detected by many 
experiments (e.g., Kovac et al. 2002; Leitch et al. 2005; 
Montroy et al. 2006; Sievers et al. 2007; Wu et al. 2007; 
Bischoff et al. 2008; Larson et al. 2010). These mea- 
surements are consistent with predictions from CMB- 
temperature-anisotropy measurements, and they provide 
new information on the epoch of reionization. Only BI- 
CEP has accurately measured E-mode polarization in 
the region of the first acoustic peak (Chiang et al. 2010); 
that paper also reports the best limit on r coming from 
cosmological B modes: r < 0.72 at the 95% confidence 
level. 

Experiments measuring B-mode polarization in the 
CMB should yield the best information on r, but this 
technique is still in its infancy. B modes are expected 
to be at least an order of magnitude smaller than the 
E modes so control of systematic errors and foregrounds 
will be particularly critical. Below « 90 GHz, the domi- 
nant foreground comes from Galactic synchrotron emis- 
sion, while at higher frequencies, emission from thermal 
dust dominates. Most planned or operating CMB polar- 
ization experiments employ bolometric detectors observ- 
ing most comfortably at frequencies > 90 GHz, so they 
cannot estimate synchrotron contamination from their 
own data. 
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26 Current address: Department of Physics, McGill University, 
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The Q/U Imaging ExpcrimcnT (QUIET) is one of two 
CMB polarization experiments to observe at frequencies 
suitable for addressing synchrotron contamination, mak- 
ing observations at 43 GHz (Q band) and 95 GHz (W 
band) and with sufficient sensitivity to begin to probe 
primordial B modes. The other is Planck (Tauber et al. 
2010). 

QUIET uses compact polarization-sensitive modules 
based upon High-Electron-Mobility Transistor (HEMT) 
amplifiers, combined with a new time-stream "double- 
demodulation" technique, side-fed Dragoman optics (for 
the first time in a CMB polarization experiment), natu- 
ral sky rotation, and frequent rotation about the optical 
axis to achieve a very low level of contamination in the 
multipolc range where a primordial-B-mode signal is ex- 
pected. 

Between 2008 October and 2010 December, QUIET 
collected over 10,000 hours of data, split between the Q- 
band and W-band receivers. Here we report first results 
from the first season of 3458 hours of Q-band observa- 
tion. The principal investigator for QUIET was our re- 
cently deceased colleague, Bruce Winstein, whose intel- 
lectual and scientific guidance were crucial to QUIET 
in all its stages, from design through analysis, through 
writing this paper. 

After describing the instrument, observations, and de- 
tector calibrations (Sections 2, 3, and 4), we discuss our 
analysis techniques and consistency checks (5 and 6). 
CMB power spectra are then presented together with 
a foreground detection (7). We evaluate our systematic 
errors (8) and then conclude (9). 

2. THE INSTRUMENT 

The QUIET instrument comprises an array of correla- 
tion polarimeters cooled to 20 K and coupled to a dual- 
rcflcctor telescope, installed on a three-axis mount in- 
side a comoving ground screen. The instrument is illus- 
trated in Figure 1. Further details are given below and in 
Newburgh et al. (2010), Kusaka et al. (2010), and Buder 
(2010). 




FlG. 1. — Overview of the QUIET instrument. The cryostat and 
1.4-m telescope mirrors are enclosed in a rectangular comoving ab- 
sorbing ground screen; in this figure its walls are transparent. The 
telescope, cryostat and electronics are mounted on a single plat- 
form attached to the deck bearing, which allows rotations around 
the instrument's optical axis. 

The Q-band QUIET receiver is a 19-element array con- 
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taining 17 low-noise correlation polarimcters, each simul- 
taneously measuring the Stokes Q, U, and I parameters, 
and two CMB differential-temperature monitors. 

QUIET uses a 1.4-m classical side- fed Dragonian an- 
tenna (Dragone 1978). This consists of a parabolic 
primary, a concave hyperbolic secondary along with a 
platelet array of corrugated feed horns (Gundersen & 
Wollack 2009). These elements are oriented in a way to 
satisfy the Mizuguchi condition (Mizugutch et al. 1976) 
in order to minimize cross polar response, and unlike dual 
offset classical Gregorian or Cassegrain antennas, the el- 
ements combine to generate high gain with low sidclobc 
response over a wide field of view (Chang & Prata 2004). 
The telescope is described in detail in Imbriale et al. 
(2010). It yields a full- width half-maximum (FWHM) 
beam size of 27f3 and a roughly circular field of view 
of 7° diameter. Radiation from each feed horn enters a 
septum polarizer (Borncmann & Labay 1995) which sep- 
arates left and right circularly-polarized components (L 
and R) into two waveguide ports which mate to a QUIET 
correlation module, detailed below. 

The module array and feed horns are cooled to 20 K in 
a cryostat to reduce instrumental noise. An electronics 
enclosure mounted next to the cryostat houses the elec- 
tronics necessary for biasing the modules and recording 
their data. The cryostat, electronics, and telescope are 
installed on the former CBI mount (Padin et al. 2002). 
This mount provides three-axis motion: azimuth, eleva- 
tion, and rotation about the optical axis. This last is 
called "deck" rotation. 

The cryostat and telescope are enclosed by an absorb- 
ing comoving ground screen. The ground screen was 
designed to have two parts, but the upper section (not 
shown in Fig. 1) was not installed until after the Q-band 
instrument was removed. Its absence was correctly antic- 
ipated to result in two far sidelobes, which were mapped 
with a high-power source by the QUIET W-band instru- 
ment in the field and measured to be < — 60 dB with 
the QUIET Q-band instrument when the Sun passed 
through them. The effects of these sidelobes are miti- 
gated through filtering and data selection (Sections 5.1.3 
and 5.2). Section 8.4 shows that any residual contami- 
nation is small. 

Each QUIET Q-band correlation module, in a foot- 
print of only 5.1 x 5.1cm 2 , receives the circular po- 
larization modes of the celestial radiation and outputs 
Stokes Q, U and / as follows. Each input is indepen- 
dently amplified and passed through a phase switch. 
One phase switch alternates the sign of the signal volt- 
age at 4 kHz, while the other switches at 50 Hz. The 
two signals are combined in a 180° hybrid coupler, with 
outputs proportional to the sum and difference of the 
inputs. Since the module inputs are proportional to 
(L, R) = (E x ±iE y ) /y/2, where E x and E y are orthogonal 
components of the incident electric field, the coupler out- 
puts are amplified versions of E x and iE y , with the phase 
switch reversing their roles. Half of each output is band- 
pass filtered and rectified by a pair of detector diodes, 
while the other half passes into a 90° hybrid coupler. A 
second pair of bandpass filters and detector diodes mea- 
sures the power from this coupler's outputs (Kangaslahti 
et al. 2006). 

Synchronous demodulation of the 4-kHz phase switch- 



ing yields measurements of Stokes +Q and — Q on the 
first two diodes and Stokes +U and — U on the remain- 
ing two. This high-frequency differencing suppresses low- 
frequency atmospheric fluctuations as well as 1/f noise 
from the amplifiers, detector diodes, bias electronics, and 
data-acquisition electronics. Subsequent demodulation 
of the 50-Hz phase switching removes spurious instru- 
mental polarization generated by unequal transmission 
coefficients in the phase-switch circuits. The resulting 
four "double-demodulated" time streams are the polar- 
ization channels. 

Averaging the output of each diode rather than demod- 
ulating it results in a measurement of Stokes /, hereafter 
called total power, denoted "TP." The TP time streams 
are useful for monitoring the weather and the stability 
of the detector responsivities, but suffer too much con- 
tamination from 1/f noise to constrain the CMB tem- 
perature anisotropy. Therefore, the Q-band instrument 
includes two correlation modules that are coupled to a 
pair of neighboring feed horns to measure the temper- 
ature difference between them, in a scheme similar to 
the WMAP differencing assemblies (Jarosik et al. 2003). 
These differential-temperature modules provide calibra- 
tion data for the telescope pointing, beams, and side- 
lobes, as well as CMB data. Their feed horns are in the 
outer ring of the close-packed hexagonal array, 3° from 
the center. 

Here we summarize several array-wide characteristics 
of the polarimcters. Bandpass measurements in the lab 
and at the start of the observing season find that the 
average center frequency is 43.1 ±0.4 GHz, and the aver- 
age bandwidth is 7.6 ± 0.5 GHz. We calculate the noise 
power spectra of the double-demodulated polarimeter 
time streams from each 40-90-minute observation to as- 
sess their 1/f knee frequencies and white-noise levels (see 
Section 5.1). The median 1/f knee frequency is 5.5 mHz, 
well below the telescope scan frequencies of 45-100 mHz. 

From the white-noise levels and responsivities (Sec- 
tion 4.1) we find an array sensitivity 28 to CMB fluctua- 
tions of 69 /iKy/s, such that the mean polarized sensitiv- 
ity per module is 280 ^Ka/s. 




Fig. 2. — The CMB and Galactic patches, in equatorial coordi- 
nates, superimposed on a Q-band all-sky WMAP 7-year tempera- 
ture map (Jarosik et al. 2010). Note that the Galactic-plane tem- 
perature signal saturates the color scale. Patch G-2 is the Galactic 
center. 



28 This is the sensitivity for 62 polarization channels. Six of 68 
polarization channels are non-functional — an array yield of 92%. 
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3. OBSERVATIONS 

QUIET is located on the Chajnantor plateau in 
the Atacama Desert of northern Chile (67°45'42"W, 
23°01'42"S). A combination of high altitude (5080m) 
and extreme dryness results in excellent observing con- 
ditions for most of the year. During the eight months 
of QUIET Q-band observations, the median prccipitablc 
water vapor (PWV) measured at the nearby APEX site 
(Glisten et al. 2006) was 1.2 mm. 

We began observations with the Q-band receiver on 

2008 October 24, and took 3458 hours of data until 

2009 June 13 (when the receiver was replaced on the 
telescope by the 90-element W-band receiver). Of these 
data, 77% are for CMB, with 12% of the observing time 
used for Galactic fields, 7% for calibration sources, and 
4% cut due to obvious instrumental problems such as lack 
of telescope motion. We observe 24 hours a day, except 
when interrupted. Our full-season operating efficiency is 
63%; causes of downtime include occasional snow, power 
outages, and mechanical failures. 



TABLE 1 

Patch Locations and Integration Times 



Patch 


RA 


Dec. 


Integration 




(J2000) 


Hours 


CMB-1 


12' l 04 m 


-39°00' 


905 


CMB-2 


05 h 12 m 


-39°00' 


703 


CMB-3 


00 h 48 m 


-48°00' 


837 


CMB-4 


22 h 44 m 


-36°00' 


223 


G-l 


16 ft 00 m 


-53°00' 


311 


G-2 


17 h 46 m 


-28° 56' 


92 



Note. — The central equatorial coordinates and integration 
times for each observing patch. G-l and G-2 are Galactic patches. 

3.1. Field Selection 

We observe four CMB fields, referred to henceforth as 
"patches." Table 1 lists their center positions and to- 
tal integration times, while Figure 2 indicates their posi- 
tions on the sky 29 . The number of patches is determined 
by the requirement to always have one patch above the 
lower elevation limit of the mount (43°). The specific 
positions of each patch were chosen to minimize fore- 
ground emission using WMAP 3-year data. The area of 
each patch is « 250 deg 2 . In addition to the four CMB 
patches, we observe two Galactic patches. These allow 
us to constrain the spectral properties of the polarized 
low-frequency foregrounds with a high signal-to-noise ra- 
tio. The results from the Galactic observations will be 
presented in a future publication. 

3.2. Observing Strategy 

Scanning the telescope modulates the signal from the 
sky, converting CMB angular scales into frequencies in 
the polarimctcr time streams. Since QUIET targets large 
angular scales, fast scanning (ss 5° s _1 in azimuth) is crit- 
ical to ensuring that the polarization modes of interest 
appear at higher frequencies than the atmospheric and 
instrumental 1// knee frequencies. 

29 Patch CMB-3 partially overlaps with the field the BICEP 
collaboration has observed for CMB analysis (Chiang et al. 2010). 
The data may be used for future analysis cross-correlating maps 
from the two experiments. 



So that each module sees a roughly-constant atmo- 
spheric signal, QUIET scans are periodic motions solely 
in azimuth with both the elevation and deck-rotation 
axes fixed. Each scan has an amplitude of 7.5° on the sky, 
with period 10-22 s. These azimuth scans are repeated 
for 40-90 minutes; each set of scans at fixed elevation is 
denoted a "constant-elevation scan" (CES). We repoint 
the telescope and begin a new CES when the patch cen- 
ter has moved by 15° in order to build up data over an 
area of « 15° x 15° for each patch. Note that a cen- 
tral region ~ 8° across is observed by all polarimeters 
since the instrument's field of view has a diameter of 
~ 7°. Diurnal sky rotation and weekly deck rotations 
provide uniform parallactic-angle coverage of the patch, 
and ensure that its peripheral regions are also observed 
by multiple polarimeters. 



TABLE 2 
Regular Calibration Observations 



Source 


Schedule 


Duration (min.) 


sky dips 


every 1.5 hours 


3 


Tau A 


every 1-2 days 


20 


Moon 


weekly 


60 


Jupiter 


weekly 


20 


Venus 


weekly 


20 


RCW38 


weekly 


20 



4. CALIBRATION 

Four quantities are required to convert polarimeter 
time streams into polarization power spectra: detector 
rcsponsivitics, a pointing model, detector polarization 
angles, and beam profiles. To this end, a suite of cali- 
bration observations is performed throughout the season 
using astronomical sources (Taurus A-hcreafter Tau A, 
Jupiter, Venus, RCW38, and the Moon); atmospheric 
measurements ("sky dips," which typically consist of 
three elevation nods of ±3°); and instrumental sources 
(a rotating sparse wire grid and a polarized broadband 
noise source). From these we also measure instrumen- 
tal polarization, as described below. QUIET's regular 
calibration observations are summarized in Table 2. 

We typically use two or more methods to determine a 
calibration constant, taking the spread among the meth- 
ods as an indication of the uncertainty. We show in Sec- 
tion 8 that aside from the case of absolute responsivity, 
all calibration uncertainties lead to estimates of system- 
atic effects on the power spectra well below statistical er- 
rors. This immunity comes from having a large number 
of detectors and highly-crosslinked polarization maps. 

4.1. Responsivity 

The polarized flux from Tau A provides a 5 mK signal 
which we observe at four parallactic angles. The sinu- 
soidal modulation of the signal induced by the changing 
parallactic angles is fitted to yield responsivity coeffi- 
cients for each detector. Figure 3 shows the response of 
the four polarization channels from the central feed horn 
to Tau A. A typical responsivity is 2.3 mV K _1 , with a 
precision from a single set of observations of 6%. The 
absolute responsivity from Tau A was measured most 
frequently for the central feed horn. We choose its +Q 
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Fig. 3. — Polarimcter responses from the central feed horn to 
the polarization of Tau A at four deck angles. The horizontal axis 
corresponds to the position angle of the receiver focal plane in 
equatorial coordinates. These data were collected with one corre- 
lation module in about 20 minutes. The errors are smaller than the 
points. From top to bottom, responses are shown for the detector 
diodes sensitive to the Stokes parameters +Q, — Q, +U, and — U, 
respectively. For each, the fitted model is plotted as a dashed line. 

diode detector to provide the fiducial absolute responsiv- 
ity. 

The responsivities of other detectors relative to the 
fiducial detector are determined with the sky dips as 
described below. We have three independent means of 
assessing the relative responsivities among polarimeters: 
from nearly-simultaneous measurements of the Moon, 
from simultaneous measurements of responses to the ro- 
tating sparse wire grid in post-season tests, and from 
Tau A measurements. The errors from these methods 
are 4%, 2%, and 6% respectively, while the error from 
the sky-dip method is 4%. All the methods agree within 
errors. 

Sky dips generate temperature signals of several 
100 mK and thus permit measurement of the TP re- 
sponsivities. The signals vary slightly with PWV. We 
estimate the slope from the data as 4% mm -1 and cor- 
rect for it. This slope is consistent with the atmospheric 
model of Pardo et al. (2001). Because the ratios of the 
responsivities for the TP and polarized signals from each 
detector diode are stable quantities within a few per- 
cent of unity, we use sky dips performed immediately 
before each CES to correct short-term variations in the 
polarimeter responsivities. The responsivities vary by 
< 10% over the course of a day, due to changing thermal 
conditions for the bias electronics. Further post-season 
tests provide a physical model: the relevant tempera- 
tures are varied intentionally while the responsivities are 
measured with sky dips. We confirm the results with the 
polarized broadband source. 

We bound the uncertainty in the absolute responsivity 
of the polarimeter array at 6%. The largest contribu- 
tions to this estimate are uncertainties in (1) the beam 
solid angle (4%, see below), (2) the response difference 
between polarized and TP signals for each diode detec- 
tor (3%), and (3) the Tau A flux (3%, Weiland et al. 
2010). The first enters in converting the flux of Tau A 
into /nK, while the second enters because although one 
fiducial diode detector is calibrated directly from Tau A, 



for the rest we find relative responsivities from sky dips 
and normalize by the fiducial diode's responsivity. 

For the differential-temperature modules, all detectors 
observe the signal from Jupiter simultaneously, providing 
the absolute responsivity for all channels upon compar- 
ison with the Jupiter flux from Weiland et al. (2010). 
Observations of Venus (Hafez et al. 2008) and RCW38 
agree with the Jupiter measurements within errors, and 
sky dips track short-term variations. We calibrate the 
absolute responsivity with 5% accuracy. 
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Fig. 4. — Map of the polarization of the Moon from one detector 
diode. The amplitude of the quadrupole polarization visible here 
is 400 mK. Similar maps are produced for all 17 polarization 
modules in the array with a single hour-long observation. The 
dotted line indicates the polarization orientation of the detector. 
Contours are spaced at intervals of 100 mK, with negative contours 
indicated by dashed lines. 



4.2. Pointing 

The global pointing solution derives from a physical 
model of the 3-axis mount and telescope tied to obser- 
vations of the Moon with the central feed horn in the 
array, as well as Jupiter and Venus with the differential- 
temperature feed horns. Optical observations are taken 
regularly with a co-aligned star camera and used to mon- 
itor the time evolution of the pointing model. 

During the first two months in the season, a mechanical 
problem with the deck-angle encoder results in pointing 
shifts. The problem was subsequently repaired. Based 
on pointing observations of the Moon and other astro- 
nomical sources, we verify that these encoder shifts are 
less than 2°. Systematic uncertainties induced by this 
problem are discussed in Section 8.1. 

After the deck-angle problem is fixed, no significant 
evolution of the pointing model is found. The differ- 
ence in the mean pointing solution between the start and 
the end of the season is smaller than 1'. Observations of 
the Moon and Jupiter also provide the relative pointing 
among the feed horns. The root mean square (RMS) 
pointing error in the maps is 3'5. 

4.3. Detector Polarization Angles 

Our primary measurement of the polarization angle 
for each detector comes from observing the radial po- 
larization of the Moon, as illustrated in Figure 4. The 
measurement has a high signal-to-noise ratio and its in- 
accuracy is dominated by systematic error due to the 
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temperature gradient of the Moon surface. One can see 
the effect in the different amplitudes of the two positive 
envelopes in Figure 4. The fluctuations of the detector 
polarization-angle measurements over many observations 
with different phases of the Moon and telescope orien- 
tations are typically 1° in RMS. Although simulations 
suggest these fluctuations can be due to the failure to ac- 
count in analysis for the temperature gradient, we conser- 
vatively assign them as upper limits on the fluctuations 
of the polarization angles during the season. Based on 
this conservative limit, we estimate the systematic error 
in the CMB-power-spectra measurement in Section 8.2, 
resulting in a negligible contribution. 

Two other less precise methods also give estimates of 
the detector angles: fits to the Tan A data, and determi- 
nation of the phases of the sinusoidal responses of all the 
detectors to rotation of the sparse wire grid. In each case, 
the differences between the detector angles determined 
by the secondary method and the Moon are described by 
a standard deviation of « 3°. However, we find a mean 
shift between the Tau A-derived and Moon-derived an- 
gles of 1?7. To estimate the errors in the angles in light 
of this shift, we use an empirical approach: in Section 8.2 
we estimate the impact on the power spectra from using 
the Tau A results instead of the Moon results, and find 
it to be small. 




I 

Fig. 5. — Top panel: Polarization beam profile from Tau A 
observations with the central feed horn. The data are overplottcd 
with the expansion in Gauss-Hcrmite polynomials described in the 
text. Bottom panel: Beam window function with errors shown by 
the gray band. 



4.4. Beam Profile and Window Function 

The polarization and differential-temperature beams 
are obtained from maps created using the full data sets of 
Tau A and Jupiter observations respectively, with square 
pixels of 1 '8 on a side. For polarization, this process 
produces the main and leakage beam maps simultane- 
ously, with the latter describing the instrumental po- 



larization. The average FWHM for the beams across 
the array is 27'3, measured with O'l precision for the 
central feed horn and for the differential-temperature 
feed horns at the edge of the focal plane. The non- 
central-polarization-horn FWHMs are measured less fre- 
quently and thus are less precisely known, with an un- 
certainty of 1'5. The beam elongation is typically small 
(1%), and its effect is further reduced by the diurnal 
sky rotation and weekly deck rotations which result in 
a symmetrized effective beam in the CMB maps. We 
compute 1-dimcnsional symmetrized beam profiles, with 
a resolution of 0'6. These profiles are modeled as a 
sum of six even Gauss-Hermite terms (Monsalve 2010). 
The main-beam solid angles are computed by integrat- 
ing these models out to 54' (roughly — 28 dB), yielding 
78.0 ± 0.4 yitsr for the differential-temperature horns and 
74.3 ± 0.7 ^sr for the central horn. An average gives 
76 [isv for all horns in the array. We also examine alter- 
native estimates such as integrating the raw beam map 
instead of the analytical fit. We assign a systematic un- 
certainty of 4% based on the differences among these 
different estimates. The systematic error includes possi- 
ble contributions from sidelobes, which we constrain to 
0.7 ± 0.4 [isv with antenna range measurements carried 
out before the observation season. 

The window functions, encoding the effect of the finite 
resolution of the instrument on the power spectra, are 
computed from the central-horn and the temperature- 
horn-profile models. The central-horn beam profile and 
window function are shown in Figure 5. The uncer- 
tainty accounts for statistical error and differences be- 
tween polarization and differential-temperature beams, 
as described in Section 8.1. 

4.5. Instrumental Polarization 

Instrumental imperfections can lead to a spurious po- 
larization signal proportional to the unpolarized CMB 
temperature anisotropy. We call this the I to Q (or U) 
leakage term. In our instrument, a fraction of the power 
input on one port of the correlation module is reflected 
because of a bandpass mismatch to the septum polar- 
izer, and a fraction of the reflected power re-enters the 
other port. The dominant monopolc term comes from 
this effect. We measure the monopolc term from the po- 
larimeter responses to temperature changes, using sky 
dips; Moon, Tau A, and Galactic signals; as well as vari- 
ations from the weather. The average magnitude is 1.0% 
(0.2%) for the Q (U) diodes. Note that the discrepancy in 
the Q and U averages was predicted from measurements 
of the properties of the septum polarizers and confirmed 
in the field. We do not correct for this effect but assign 
systematic errors as described in Section 8.3. 

5. DATA ANALYSIS PROCEDURE 

QUIET employs two independent analysis pipelines to 
derive CMB power spectra. We present the methods used 
for analysis in each pipeline, including data selection, 
filtering, map making, and power-spectra estimation. 

Pipeline A is based on the pseudo-C^ analysis frame- 
work, first described by Hivon et al. (2002), which is used 
by numerous experiments (Netterfield et al. 2002; Brown 
et al. 2009; Chiang et al. 2010; Larson et al. 2010; Lueker 
et al. 2010). This pipeline made all analysis choices in 
accordance with a strict (blind) analysis validation policy 
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described in Section 6. An advantage of the pseudo-Cf 
framework is computational efficiency, which is critical 
for completing the more than 30 iterations of the null- 
test suite. For the same reason, this pipeline is used 
for the systematic-error evaluations found in Section 8. 
Pscudo-C^ analysis also enables us to perform cross cor- 
relation, making the resultant power spectra immune to 
possible misestimation of noise bias. 

Pipeline B implements a maximum-likelihood frame- 
work (e.g., Tegmark 1997; Bond et al. 1998), which has a 
long history of use by CMB experiments (e.g., Mauskopf 
et al. 2000; Page ct al. 2007; Wu ct al. 2007; Bischoff 
et al. 2008). This framework yields minimum- variance 
estimates of the power spectra, naturally accounts for 
E/B mixing, and directly provides the exact CMB like- 
lihood required for estimation of cosmological parame- 
ters, without the use of analytical approximations. In 
addition to power spectra, it produces unbiased maps 
with full noise-covariance matrices, useful for compar- 
isons with other experiments. On the other hand, this 
approach is also computationally more expensive than 
the pseudo-C^ framework, and a reduced set of null tests 
is therefore used to evaluate data consistency. 

The processing of the time-ordered data (TOD) and 
the methodology used for data selection are treated in 
Sections 5.1 and 5.2, respectively Brief descriptions of 
the pseudo-C^ and maximum-likelihood techniques are 
found in Section 5.3. TOD processing, data selection, 
and analysis for temperature-sensitive modules are dis- 
cussed in Section 5.4. 



5.1. Time- Ordered- Data Processing 

To prepare the TOD for map making, we execute three 
steps: pre-processing, noise modeling, and filtering. Of 
these steps, only the filtering is significantly different be- 
tween the two pipelines. 

5.1.1. Pre-processing 

The first data-processing step is to correct for a small 
non-linearity that was discovered in the analog-to-digital 
converter (ADC) system. The non-linearities occur ev- 
ery 1024 bits; roughly 14% of the data are affected. Sys- 
tematic uncertainty from this effect is estimated in Sec- 
tion 8.5. Next, the receiver data are synchronized with 
the telescope pointing. The double-demodulation step, 
described in Section 2, is applied, reducing the sample 
rate from 100 Hz to 50 Hz. A model of the detectors' po- 
larized responsivities converts the data from ADC counts 
into thermodynamic temperature. The two pipelines use 
different responsivity models. Pipeline A applies a con- 
stant rcsponsivity throughout each CES, addressing pos- 
sible variability within a CES as part of the systematic 
error (Section 8); pipeline B updates responsivities on 
2-minute timescales (Dumoulin 2010). 

5.1.2. Noise Model 

After pre-processing, the time streams for each detec- 
tor diode in each CES arc Fourier-transformed and their 



noise power spectra are fit to a model 30 with three pa- 
rameters: the amplitude of white noise, the 1// knee 
frequency, and the power-law slope of the 1// noise. We 
also compute the white-noise correlations among detec- 
tor diodes in the same module: the most important are 
between the two Q or the two U detector diodes (with 
an average coefficient of 0.22). A small fraction of the 
noise spectra contain features not accounted for in the 
noise model: beam sidclobes (see Section 2) scanning 
across features on the ground create a narrow spike at the 
scan frequency; slowly-changing weather patterns during 
a CES create a broader peak also at the scan frequency; 
and there are some narrow spikes at high (> 6 Hz) fre- 
quencies. To prevent these features from biasing the 
noise model, the fit excludes a region around the scan 
frequency as well as frequencies above 4.6 Hz. In ad- 
dition to the noise-model parameters, several statistics 
quantifying the agreement between the data and noise 
model are also used for data selection as described in 
Section 5.2. 

5.1.3. Filtering 

In pipeline A, three filters are applied. These were cho- 
sen from the results of many runs of the null-test suite 
(see Section 6). First, to remove the high-frequency nar- 
row spikes, we apply a low-pass filter that cuts signals 
off sharply above 4.6 Hz 31 . Second, to suppress contam- 
ination from atmospheric fluctuations and detector 1/f 
noise, we subtract a linear function from each telescope 
half scan (left-going or right-going) removing modes be- 
low twice the scan frequency 32 . The third filter, designed 
to eliminate signal from ground emission, removes any 
azimuthal structure that remains after summing over all 
half scans in the CES. 

In pipeline B, an apodized bandpass filter is used that 
accepts modes from 2.5 times the scan frequency to 
4.5 Hz; the highpass component of this filter is designed 
to suppress scan-synchronous contamination. Further, a 
time-independent ground-emission model is subtracted. 
The model of ground emission is generated by building 
low-resolution and high-signal-to-noise maps in horizon 
coordinates from the full-season data for each deck angle 
and module, using large (55') pixels. Only features that 
are stable in time, azimuth, elevation, and deck angle 
contribute to this model. The amplitude of the ground 
correction is < 1 /iK. 

5.2. Data Selection 

The fundamental unit of data used for analysis is the 
double-demodulated output of one detector diode for a 
single CES, referred to as a "CES-diode." Selecting only 
those CES-diodes that correspond to good detector per- 
formance and observing conditions is a critical aspect of 
the data analysis. The data-selection criteria began with 
a nominal set of cuts and evolved into several distinct 
configurations, as many as 33 in the case of pipeline A. 
For each configuration, analysis validation (see Section 6) 
was performed yielding statistics quantifying the lack of 

30 At the level of a single CES, the TOD of each detector diode 
are dominated by noise; the contribution of the CMB is negligible. 

31 For QUIET's beam size and scanning speed a low-pass filter 
of 4.5— 4. 6 Hz results in a minimal loss of sensitivity to the CMB. 

32 Typical scan frequencies range from 45mHz to 100 mHz. 
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contamination in the data set. The final data set was 
chosen when these statistics showed negligible contami- 
nation and were little affected by changes to the cuts. 

Cut efficiencies, defined as the fractions of CES-diodcs 
accepted for the analysis, are given for both pipelines in 
Table 3. While each pipeline applies its own cuts uni- 
formly to all four patches, the efficiencies among patches 
are non-uniform because of differences in weather qual- 
ity. Over the course of the eight month observing sea- 
son, patch CMB-1 is primarily visible at night, when the 
atmosphere tends to be more stable; patch CMB-3 is 
mostly observed during the day. 

The first step of the data selection is simply to remove 
known bad data: data from six non-functional detector 
diodes, data during periods of mount malfunctions, and 
CEScs lasting less than 1000 s. Further, we cut individual 
CES-diodes that show deviation from the expected linear 
relationship between the demodulated and TP signals. 
This cut removes data with poor thermal regulation of 
the electronics or cryostat, or residual ADC non-linearity. 

The beam sidelobes, described in Section 2, introduce 
contamination to the data if the telescope scanning mo- 
tion causes them to pass over the ground or the Sun. 
Ground pickup is dealt with by filtering as described in 
Section 5.1.3. The less frequent cases of Sun contamina- 
tion are handled by cutting those CES-diodes for which 
the Sun's position overlaps with the measured sidelobc 
regions for each diode. 

Additional cuts are specific to each pipeline. Pipeline 
A removes data taken during bad weather using a statis- 
tic calculated from fluctuations of the TP data during 
10-s periods, averaged across the array. This cut removes 
entire CESes. Several more cuts remove individual CES- 
diodes. While these additional cuts are derived from the 
noise modeling statistics, they also target residual bad 
weather. During such marginal weather conditions only 
some channels need to be cut, since the sensitivity for 
a given detector diode to atmospheric fluctuations de- 
pends on its level of instrumental polarization. Next, we 
reject CES-diodes with poor agreement between the fil- 
tered data and the noise model in three frequency ranges: 
a narrow range (only 40 Fourier modes) about the scan 
frequency, from twice the scan frequency to 1Hz, and 
from 1 Hz to 4.6 Hz. We also cut CES-diodes that have 
higher than usual 1/f knee frequencies, or large vari- 
ations during the CES in the azimuthal slopes of the 
double-demodulated time streams; both these cuts help 
eliminate bad weather periods. Finally, we also remove 
any CES-diodcs with an outlier greater than 6 a in the 
time domain on three timescales (20 ms, 100 ms, and 1 s). 

For pipeline B, the weather cut rejects CESes based 
on a statistic computed from fluctuations of the double- 
demodulated signals from the polarization modules on 
10-s and 30-s timescales. Three cuts are applied to re- 
move individual CES-diodes. The first is a cut on the 1 // 
knee frequency, similar to that of pipeline A. Second, a 
cut is made on the noise model x 2 m the frequency range 
passed by the filter, and third, we reject CES-diodes hav- 
ing a large x 2 m the azimuth-binned TOD. This cut re- 
jects data with possible time variation in the ground sig- 
nal. Finally, an entire CES is removed if more than 40% 
of its detectors have already been rejected. 

5.3. Map Making and Power-Spectra Estimation 



TABLE 3 

Total Hours Observed and Data-Selection Efficiencies 



Patch 


Total Hours 


A % 


B % 


Common % 


CMB-1 


905 


81.7 


84.3 


76.7 


CMB-2 


703 


67.3 


70.0 


61.2 


CMB-3 


837 


56.0 


61.4 


51.4 


CMB-4 


223 


70.6 


74.2 


65.9 


All Patches 


2668 


69.4 


72.9 


64.2 



Note. — Selection efficiencies for each pipeline. "Common" 
gives the efficiencies if both sets of cuts were applied. 

After filtering, the TOD for all diodes are combined to 
produce Q and U maps for each of the QUIET patches. 
The maps use a HEALPix N s jd c = 256 pixclization 
(Gorski et al. 2005). This section describes the map mak- 
ing and power-spectra estimation from the maps for each 
of the pipelines. 

5.3.1. Pipeline- A Map Making 

Polarization maps (Q and U) are made by summing 
samples into each pixel weighted by their inverse vari- 
ance, calculated from the white-noise amplitudes. The 
full covariance matrix is not calculated. Two polarized 
sources, Centaurus A and Pictor A, are visible in the 
maps and are removed using circular top-hat masks with 
radii of 2° and 1°, respectively. 

Separate maps are made for each range of telescope 
azimuth and deck-angle orientations. The coordinates 
are binned such that there are 10 divisions in azimuth 33 
and six distinct ranges of deck-angle orientation. Making 
separate maps for different telescope pointings enables 
the cross correlation described in the next section. 

5.3.2. Power-Spectra Estimation in Pipeline A 

The MASTER (Monte Carlo Apodizcd Spherical 
Transform Estimator) method is used in pipeline A 
(Hivon et al. 2002; Hansen & Gorski 2003); it is based 
on a pseudo-C^ technique and takes account of effects 
induced by the data processing using Monte Carlo (MC) 
simulations. The pscudo-C^ method allows estimation 
of the underlying Ce using spherical-harmonics transfor- 
mations when the observations do not cover the full sky 
uniformly (Wandclt et al. 2001). The pscudo-C^ spec- 
trum, designated by Ci, is related to the true spectrum 
C t by: 

(C t ) = Y,M tt >F t ,Bl{C e ). (1) 

v 

There is no term corresponding to noise bias, which 
would arise if we did not employ a cross-correlation tech- 
nique. Here E>£ is the beam window function, described 
in Section 4.4, and Mw is a mode-mode-coupling kernel 
describing the effect of observing only a small fraction 
of the sky with non-uniform coverage. It is calculable 
from the pixel weights, which are chosen to maximize 
the signal-to- noise ratio (Feldman et al. 1994). We bin 
in I and recover Ci in nine band powers, Cb, and Ft is the 
transfer function (displayed in Section 7) due to filtering 
of the data; its binned estimate, Ff,, is found by process- 
ing noiseless CMB simulations through pipeline A and 

33 The azimuth divisions are the same for all patches, which 
means that not all divisions arc populated for patches CMB-3 and 
CMB-4. 
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used to obtain Cb- For the polarization power spectra, 
equation (1) is generalized for the case where Ce contains 
both Cf E and Cf B . 

In the power-spectra estimates, we include only the 
cross correlations among pointing-division maps, exclud- 
ing the auto correlations. Because the noise is uncorre- 
lated for different pointing divisions, the cross-correlation 
technique allows us to eliminate the noise-bias term and 
thus the possible residual bias due to its misestimate. 
Cross correlation between different pointing divisions 
also suppresses possible effects of ground contamination 
and/or time- varying effects. Dropping the auto correla- 
tions creates only a small increase in the statistical errors 
(« 3%) on the power spectra. 

The errors estimated for the pipeline-A power spec- 
tra are frequentist two-sided 68% confidence intervals. A 
likelihood function used to compute the confidence in- 
tervals is modeled following Hamimeche & Lewis (2008) 
and calibrated using the MC simulation ensemble of more 
than 2000 realizations with and without CMB signal. We 
also use the likelihood function to put constraints on r 
and calculate the consistency to ACDM. 

The partial sky coverage of QUIET generates a small 
amount of E/B mixing (Challinor & Chon 2005), which 
contributes an additional variance to the BB power spec- 
trum. We incorporate it as part of the statistical error. 
This mixing can be corrected (Smith & Zaldarriaga 2007) 
in future experiments where the effect is not negligible 
compared to instrumental noise. 

5.3.3. Pipeline- B Map Making 

In pipeline B, the pixel-space sky map rh {N s a B = 256) 
is given by 



m = (P T N" 1 FP) P T N" 1 Fd, (2) 

where P is the pointing matrix, N is the TOD-noise- 
covariance matrix, F corresponds to the apodized band- 
pass filter discussed in Section 5.1.3, and d denotes the 
TOD. This map is unbiased, and for the case F = 1 it is 
additionally the maximum-likelihood map, maximizing 



£(m|d) 



-■I (d-Pm) T N" 1 (d-Pm) 



(3) 



The corresponding map-noise-covariance matrix (e.g., 
Tegmark 1997; Keskitalo et al. 2010) is 

Nrh = (P T N" 1 FP) _1 (P T F T N~ 1 FP) (P T N _1 FP) _1 . 

Note that one often encounters the simplified expression 
N A = (P^^FP) -1 in the literature. This corre- 
sponds effectively to assuming that F = F 2 in the Fourier 
domain, and is strictly valid for top-hat-filter functions 
only. For our filters, we find that the simplified expres- 
sion biases the map-domain x 2 (= f^N^n, where n is 
a noise-only map) by « 3 a, and we therefore use the full 
expression, which does lead to an unbiased \ 2 . 

Equations (2-4) apply to both polarization and tem- 
perature analysis. The only significant difference lies in 
the definition of the pointing matrix, P. For polariza- 
tion, P encodes the detector orientation, while for tem- 
perature it contains two entries per time sample, +1 and 
— 1, corresponding to the two horns in the differential- 
temperature assembly. 



After map making, the maps are post-processed by re- 
moving unwanted pixels (i.e., compact sources and low- 
signal-to- noise edge pixels). All 54 compact sources in 
the 7- year WMAP point source catalog (Gold ct al. 2010) 
present in our four patches are masked out, for a total 
of 4% of the observed area. We also marginalize over 
large-scale and unobserved modes by projecting out all 
modes with £ < 5 (£ < 25 for temperature) from the 
noise-covariance matrix using the Woodbury formula, as- 
signing infinite variance to these modes. 

5.3.4. Power- Spectra Estimation in Pipeline B 

Given the unbiased map estimate, m, and its noise- 
covariance matrix, Nj,, we estimate the binned CMB 
power spectra, Cb, using the Newton-Raphson optimiza- 
tion algorithm described by Bond et al. (1998), general- 
ized to include polarization. In this algorithm one iter- 
ates towards the maximum-likelihood spectra by means 
of a local quadratic approximation to the full likelihood. 
The iteration scheme in its simplest form is 



SC b = ^ ^Tr [(mm T - CXC^C, yC" 1 )] 



(5) 



where b denotes a multipole bin, C is the signal-plus- 
noise pixel-space covariance matrix, and is the 
derivative of C with respect to Cb- The signal compo- 
nent of C is computed from the binned power spectra, 
Cb, and the noise component is based on the noise model 
described in Section 5.1.2, including diode-diode correla- 
tions. Finally, 



Tbw = i r &(C- 1 C l6 C- 1 C 1 6') 



(6) 



is the Fisher matrix. Additionally, we introduce a step 
length multiplier, a, such that the actual step taken at 
iteration i is a 5Cb, where < a < 1 guarantees that C is 
positive definite. We adopt the diagonal elements of the 
Fisher matrix as the uncertainties on the band powers. 

We start the Newton-Raphson search at Ce — 0, and 
iterate until the change in the likelihood value is lower 
than 0.01 times the number of free parameters, corre- 
sponding roughly to a 0.01-cr uncertainty in the position 
of the multivariate peak. Typically we find that 3 to 10 
iterations arc required for convergence. 

Estimation of cosmological parameters, 9, is done by 
brute-force grid evaluation of the pixel-space likelihood, 



C{6) 



id-C 



(9)d 



(7) 



Here C(9) is the covariance matrix evaluated with a 
smooth spectrum, Ce, parametrized by 9. In this pa- 
per, we only consider 1-dimensional likelihoods with a 
parametrized spectrum of the form Ce = aCf d , a being 
a scale factor and Cf d a reference spectrum; the compu- 
tational expense is therefore not a limiting factor. Two 
different cases are considered, with a being either the 
tensor-to-scalar ratio, r, or the amplitude of the EE spec- 
trum, q, relative to the ACDM model. 

5.4. Temperature Data Selection and Analysis 

As described in Section 2, we dedicate one pair of mod- 
ules to differential-temperature measurements. While 
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these modules arc useful for calibration purposes, when 
combined with our polarization data they also enable us 
to make self-contained measurements of the TE and TB 
power spectra. 

For temperature, both pipelines adopt the pipeline- 
A data-selection criteria used for polarization analysis 
(see Section 5.2). The temperature-sensitive modules, 
however, are far more susceptible to atmospheric con- 
tamination than the polarization modules. Thus, these 
cuts result in reduced efficiencies: 12.4%, 6.9%, and 6.8% 
for patches CMB-1, CMB-2, and CMB-3, respectively 34 . 
More tailoring of the cuts for these modules would im- 
prove efficiencies. 

In pipeline A, the analysis proceeds as described in Sec- 
tions 5.1.3, 5.3.1, and 5.3.2 except for two aspects. First, 
in the TOD processing a second-order polynomial is fit 
and removed from each telescope half scan instead of a 
linear function. This suppresses the increased contami- 
nation from atmospheric fluctuations in the temperature 
data. Second, we employ an iterative map maker based 
on the algorithm described by Wright et al. (1996). Map 
making for differential receivers requires that each pixel 
is measured at multiple array pointings or crosslinkcd. 
In order to improve crosslinking wc divide the tempera- 
ture data into only four maps by azimuth and deck an- 
gle, rather than the 60 divisions used for polarization 
analysis. To calculate TE and TB power spectra, polar- 
ization maps are made for these four divisions, plus one 
additional map that contains all polarization data with 
pointings not represented in the temperature data. 

For pipeline B the algorithms for making temperature 
maps and estimating power spectra are identical to the 
polarization case, as described in Sections 5.3.3 and 5.3.4. 

6. ANALYSIS VALIDATION 

The QUIET data analysis follows a policy of not look- 
ing at the power spectra until the analysis is validated 
using a set of predefined tests for possible systematic ef- 
fects 35 . The validation tests consist of a suite of null 
tests, comparisons across multiple analysis configura- 
tions, and consistency checks among power spectra from 
different CMB patches. Data-selection criteria, filtering 
methods, and the division of data into maps for cross 
correlation in pipeline A are all evaluated based on the 
test results. 

Details of tests found in this section describe pipeline 
A. While the pipeline B analysis follows a similar pro- 
gram of null tests to verify the result, the increased 
computational requirements of the maximum-likelihood 
framework limit the number of tests that could be per- 
formed and require those tests to be run using lower- 
resolution maps than for the non-null analysis. The bulk 
of this section treats validation of the polarization analy- 
sis; at the end, we briefly describe the temperature anal- 
ysis validation. 

We conduct this validation in a blind-analysis frame- 
work to reduce experimenter bias, the influence of the 
experimenter's knowledge of prior results and theoretical 
predictions on the result (power spectra). Blind analysis, 

34 Patch CMB-4 is excluded due to low data-selection efficiency 
and a lack of sufficient crosslinking. 

35 Some systematic effects, such as a uniform responsivity- 
calibration error, cannot be detected by these techniques, and are 
addressed in Section 8. 



making the analysis choices without knowing the result, 
is a standard technique for minimizing this bias (Klein 
& Roodman 2005). In our blind analysis framework, we 
finalize all aspects of the data analysis including calibra- 
tion, data selection, and evaluation of the systematic er- 
ror. Only after the analysis is finalized and the following 
validation tests pass do we examine the power spectra. 

In a null test, the data are split into two subsets. Maps, 
mi and m2, are made from each subset. The power spec- 
tra of the difference map, m^is = (mi — m2)/2, are ana- 
lyzed for consistency with the hypothesis of zero signal. 
The null suite consists of 42 tests 36 , each targeting a 
possible source of signal contamination or miscalibration. 
These are highly independent tests; the data divisions for 
different null tests are correlated at only 8.8% on average. 
Nine tests divide the data by detector diode based on sus- 
ceptibility to instrumental effects, such as instrumental 
polarization. Ten tests target effects that depend on the 
telescope pointing such as data taken at high or low el- 
evation. Five tests divide based on the proximity of the 
main or sidelobe beams to known sources such as the Sun 
and Moon. Eight tests target residual contamination in 
the TOD using statistics mentioned in Section 5.2. Ten 
tests divide the data by environmental conditions such 
as ambient temperature or humidity. 

Each null test yields EE and BB power spectra in nine 
i bins, calculated separately for each CMB patch. Figure 
6 shows the power spectra from one null test. Although 
the EB spectra are also calculated for each null test, they 
are assigned lesser significance since sources of spurious 
EB power will also result in the failure of EE and BB null 
tests. Combining all EE and BB points for all patches 
and null tests in the null suite yields a total of 3006 null- 
spectrum points. For each power-spectrum bin 6, wc 
calculate the statistic Xnuii = C™ n /cr(,, where C™ 11 is the 
null power and Ob is the standard deviation of C™ 11 in 
MC simulations. We evaluate both Xnuii and its square; 
Xnuii is sensitive to systematic biases in the null spectra 
while Xnuii is more responsive to outliers. We run MC 
simulations of the full null suite to take into account the 
small correlation among the null tests and the slight non- 
Gaussianity of the Xnuii distribution. Non-Gaussianity is 
caused by the small number of modes at low i. 




500 



Fig. 6. — EE and BB power spectra for the patch CMB-1 null 
test between Q and U detector diodes. The inset shows the low-£ 
region in detail. 

36 Only 41 null tests are performed for patch CMB-4; one test 
is dropped because there are no data in one of the subsets. 
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As we refine the data-selection criteria based on the 
results of the null suite, we use a second test to monitor 
changes in the non-null power spectra. Using a blind 
analysis framework, we compute the difference of the 
power spectra between any two iterations of the data se- 
lection without revealing the non-null spectra. Further, 
we randomize the sign of the difference to hide the direc- 
tion of the change; knowledge of the direction could allow 
experimenter bias (e.g. a preference for low BB power). 
Figure 7 shows the differences in the power spectra be- 
tween the final configuration and several intermediate 
iterations of the data selection, starting with data sets 
that showed significant failures for the null-test suite. 
Statistically significant differences indicate a change in 
the level of contamination in the selected data set. Our 
data-selection criteria are finalized when further itera- 
tions only result in statistically expected fluctuations. 
The sensitivity of this test is demonstrated by the fact 
that the expected fluctuations are much less than the 
statistical error of the final result. 

Finally, the non-null power spectra are compared 
among the four CMB patches. A % 2 statistic is computed 
from the deviation of each patch's non-null power spectra 
from the weighted average over all patches. The total x 2 
is compared to MC simulations to compute probabilities 
to exceed (PTE). 
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Fig. 7. — Power-spectra differences between the final data se- 
lection and six of the 32 earlier data-selection iterations, ordered 
by date. The lowest-^ bin of patch CMB-1 is shown. The error 
bars correspond to the expected fluctuations due to the differences 
in data selected, which are much smaller than the final statistical 
errors in this bin 0.10 /iK 2 for BB). Iterations that are closer to 
the final data selection have smaller errors. The expected EE power 
in this bin from the ACDM model is also shown for comparison. 

When all aspects of the analysis are finalized, the last 
round of null tests and CMB patch comparisons validates 
the non-null-power-spectra results. Figure 8 shows the 
distributions of the Xnuii statistic and of the PTEs cor- 
responding to all XnuXi va l ues from the full null suite. In 
pipeline A, the distribution of Xnuii is consistent with the 
expectation from MC simulations. The mean of the Xnuii 
distribution is 0.02 ±0.02; the mean of the MC-cnscmble 
Xnuii distribution is also consistent with zero. The distri- 
bution of the Xnuii PTEs is uniform as expected. Table 4 
lists the PTEs for the sums of the Xnuii statistic over all 
bins in each patch. Examinations of various subsets of 
the null suite, such as EE or BB only, do not reveal any 



anomalies. The EB null spectra do not indicate any fail- 
ure either. Patch comparison PTEs are 0.16, 0.93, and 
0.40 for EE, BB, and EB, respectively, demonstrating no 
statistically significant difference among the patches. 

A similar, but smaller, null suite is run by pipeline B. 
Specifically, 21 null tests are made at a HEALPix res- 
olution of iV s idc = 128. The results obtained in these 
calculations are summarized in the bottom panel of Fig- 
ure 8, and total PTEs for each patch are listed in Table 4. 
As in pipeline A, no anomalous values are found. 

Finally, we make a comment on the usefulness of the 
Xnuii distribution (as opposed to the Xnuii distribution) 
for identifying and quantifying potential contaminants. 
During the blind stage of the analysis, a positive bias in 
the Xnuii distribution of 0.21 (0.19) was identified using 
pipeline A (B) (corresponding to 21% (19%) of the statis- 
tical errors). The number from pipeline A was obtained 
when including auto correlations in its power-spectra es- 
timator. When excluding auto correlations, and cross- 
correlating maps made from data divided by time (day 
by day), the bias decreased to 0.10. Further detailed 
studies lead to the division of data into maps based on 
the telescope pointing, as described in Section 5.3; the 
result is an elimination of the observed bias. 

The maximum-likelihood technique employed by 
pipeline B intrinsically uses auto correlations, and a cor- 
responding shift in the Xnuii distribution is seen in Figure 
8. However, as will be seen in Section 7, the power spec- 
tra from the two pipelines are in excellent agreement, 
thereby confirming that any systematic bias coming from 
including auto correlations is well below the level of the 
statistical errors. We close this section by mentioning 
that we know of no other CMB experiment reporting an 
examination of the Xnuii distribution, which is sensitive 
to problems not detected by examining the XnuXi distri- 
bution only. 

6.1. Validation of the Temperature Analysis 

A smaller number of null tests is used for the temper- 
ature analysis. Several are not applicable and others are 
discarded due to lack of data with sufficient crosslink- 
ing. Even so, we are able to run suites of 29, 27, and 23 
TT null tests on patches CMB-1, CMB-2, and CMB-3, 
respectively. We calculate the sums of Xnuii statistics, 
yielding PTEs of 0.26 and 0.11 for patches CMB-1 and 
CMB-2, respectively. No significant outliers are found for 
these patches. However, a 5-cr outlier in a single test 37 
is found in patch CMB-3, implying contamination in its 
temperature map. CMB-3 is therefore excluded from 
further analysis. We confirm consistency between the 
patches CMB-1 and CMB-2 with a PTE of 0.26. 

With no significant contamination in TT, EE, or BB 
spectra, one may be confident that the TE and TB spec- 
tra are similarly clean. For confirmation, we calculate 
TE and TB null spectra for the five null tests that are 
common to the temperature and polarization analyses. 
These yield PTEs of 0.61 and 0.82 for TE, and 0.16 and 
0.55 for TB, for patches CMB-1 and CMB-2, respectively, 
with no significant outliers. Patch consistency checks 
give PTEs of 0.48 for TE and 0.26 for TB. Thus, the 
TE and TB power spectra, as well as the TT, pass all 
validation tests that are performed. 

37 This null test divides the data based on array pointing. 
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Fig. 8. — Null-Suite Statistics. The upper panel shows a his- 
togram of the Xnull values for the pipclinc-A null suite (circles), 
pipclinc-B null suite (triangles), and the average of 1024 MC real- 
izations of the pipclinc-A null suite (gray histogram). Both data 
and MC distributions show similar non-Gaussianity in the Xnull 
statistic. The shift in Xnull seen for pipeline B, also seen in ear- 
lier iterations of pipeline A, is discussed in the text. The lower 
panel shows a histogram of PTEs calculated from the X nu ii statis- 
tic (outliers from either side of the upper distribution manifest as 
low PTEs). 



TABLE 4 

Null Suite Probability To Exceed by Patch 



Patch 


Pipeline A % 


Pipeline B % 


CMB-1 


44 


7 


CMB-2 


19 


43 


CMB-3 


16 


23 


CMB-4 


68 


28 



Note. — PTEs calculated from the sums of the X nu ii statistics, 
for EE and BB spectra points, over the null tests for each patch. 

7. RESULTS 

We report results from the first season of QUIET Q- 
band observations: CMB power spectra, derived fore- 
ground estimates, and constraints on the tensor-to-scalar 
ratio, r. 

7.1. Polarization Power Spectra 

The CMB power spectra are reported in nine equally- 
spaced bands with A£ = 50, beginning at £ min = 25. 
Given the patch size, modes with £ < £ m i n cannot be 
measured reliably. The correlation between neighboring 
bins is typically —0.1; it becomes negligible for bins fur- 
ther apart. 

The EE, BB, and EB polarization power spectra es- 
timated by both pipelines are shown in Figure 9. The 
agreement between the results obtained by the two 
pipelines is excellent, and both are consistent with the 
ACDM concordance cosmology. Our findings and con- 
clusions are thus fully supported by both pipelines. Only 
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Fig. 9. — EE, BB, and EB power spectra from each QUIET 
pipeline, all four patches combined. The insets show the low-£ 
region in detail. Window and transfer functions for each i bin 
are shown below the corresponding power spectra in black and 
gray, respectively. The window function combines the mode-mode— 
coupling kernel M^i with the beam (Bg ) and represents, in combi- 
nation with the transfer function (Fi), the response in each band 
to the true C( spectrum. The EE point in the lowest-^ bin in- 
cludes foreground contamination from patch CMB-1. For this dis- 
play, pipeline A shows frequentist 68% confidence intervals while 
pipeline B uses the diagonal elements of the Fisher matrix; the dif- 
ference is most pronounced in the lowest-^ bin where the likelihood 
is the most non-Gaussian. 

the statistical uncertainties are shown here; we treat sys- 
tematic errors in Section 8. Because the systematic er- 
ror analysis was only done for pipeline A, we adopt its 
power-spectra results (tabulated in Table 5) as the offi- 
cial QUIET results. 

The bottom sub-panels in Figure 9 show the window 
and transfer functions for each bin computed by pipeline 
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A. Figure 10 shows the maps for patch CMB-1 com- 
puted by pipeline B, and Figure 11 shows the QUIET 
power spectra in comparison with the most relevant ex- 
periments in our multipole range 38 . Additional plots and 
data files are online 39 . 

Fitting only a free amplitude, q, to the EE spectrum 
relative to the 7-year best-fit WMAP ACDM spectrum 
(Larson et al. 2010), we find q = 0.87 ± 0.10 for pipeline 
A and q = 0.94±0.09 for pipeline B. Taking into account 
the full non-Gaussian shapes of the likelihood functions, 
both results correspond to more than a 10-cr detection of 
EE power. In particular, in the region of the first peak, 
76 < I < 175, we detect EE polarization with more than 
6-<7 significance, confirming the only other detection of 
this peak made by BICEP at higher frequencies. The % 2 
relative to the ACDM model, with Cf B = Cf B = 0, is 
31.6 (24.3) with 24 degrees of freedom, corresponding to 
a PTE of 14% (45%) for pipeline A (B). 

7.2. Foreground Analysis 

In order to minimize possible foreground contamina- 
tion, QUIET's four CMB patches were chosen to be far 
from the Galactic plane and known Galactic synchrotron 
spurs. In these regions, contributions from thermal dust 
emission are negligible in Q band. Spinning dust is ex- 
pected to be polarized at no more than a few percent 
in Q band (Battistelli et al. 2006; Lopez-Caraballo et al. 
2010), so we expect the contribution to polarized fore- 
ground emission in our patches to be small. We there- 
fore consider only two dominant sources of possible fore- 
ground contamination, namely compact radio sources 
and Galactic diffuse synchrotron emission. 

To limit the effect of compact radio sources, we apply 
a compact-source mask to our maps before computing 
the power spectra, as described in Section 5. We also 
evaluate the CMB spectra both with and without the full 
WMAP temperature compact-source mask (Gold et al. 
2010), and find no statistically significant changes. The 
possible contribution from compact radio sources with 
fluxes below the WMAP detection level (1 Jy) is small: 
0.003 ^K 2 at £ = 50 and 0.01 ^K 2 at I = 100 (Battye 
et al. 2010) 41 . We therefore conclude that our results are 
robust with respect to contamination from compact radio 
sources and that the dominant foreground contribution 
comes from diffuse synchrotron emission. 

In Figure 12 we show the power spectra measured from 
each patch. The CMB-1 EE band power for the first 
bin is 0.55 ± 0.14 /iK 2 , a 3-ct outlier relative to the ex- 
pected ACDM band power of 0.13 /xK 2 ; while not signif- 
icant enough to spoil the overall agreement among the 
patches as shown in Section 6, this is a candidate for a 
bin with foreground contamination. 

To estimate the Q-band polarized synchrotron contam- 
ination in our CMB patches, we process the WMAP1 

38 Since Larson et al. (2010) does not provide an upper limit on 
the BB power, we use the diagonal elements of the Fisher matrix 
and show the points that are 1.65 a above their central values as 
95% upper limits. 

39 http:/ /quiet. uchicago.edu/rcsults/index. html 

40 Only t > 76 are used in the EE fit and the x 2 calculation rela- 
tive to ACDM because the first EE bin has a significant foreground 
contribution; see Section 7.2. 

41 The estimate is robust since 90% of the contribution comes 
from the high-flux population between 100 mjy and 1 Jy, where the 
distribution of the population is well understood. 



TABLE 5 

CMB-Spectra Band Powers from QUIET Q-Band Data 



bin 



EE 



1313 



EB 



25-75 
76-125 
126-175 
176-225 
226-275 
276-325 
326-375 
376-425 
426-475 



J 0.33 



+0.16 
0.11 



-0.01] 



0-82±g;lg 

,+0.34 
-0.31 



0.93 

I. 11 

2.46 
8.2 

II. 5 



+0.58 
0.52 
+ 1.10 
0.99 
+2.1 
1.9 
+3.6 

+ 13 
11 



2r 



,+0.06 
-0.04 

0-04t° ;» 

0.241^ 

64+ 53 

1 n 7 +0.98 
iiU '-0.86 

+ 1.6 
1.4 
+2.7 
2.4 
+5.3 
4.9 
11 

10 



0.00 
-0.10 
0.71 



+0.07 
0.07 
+0.11 
0.12 
+0.22 
0.20 



0.8 
-2.2 
-4.9 
2 



n ,0+0.38 
U - i8 -0.38 

-o.52l° ;«« 



-,+2.0 
-2.0 
o 9 +3.9 
°- z -3.9 
;+8.3 
-8.2 



o.o 1 



4.5" 1 



Note. — Units are thermodynamic temperatures, /iK 2 , scaled 
as C e e(£ + 1)/2tt. 

a Patch CMB-1 has significant foreground contamination in the 
first EE bin. 




Stokes U 




-15 uKl 



E modes 





\ 9 






pe— ■+ 
- • 


. 




V 


280° 



-1.5nKI 




1.5uK 



Fig. 10. — Maps of patch CMB-1 in Galactic coordinates. The 
top row shows our polarization maps with compact sources masked 
(white disks). The bottom row shows E and B modes decom- 
posed using a generalized Wiener filter technique, implemented 
through Gibbs sampling where the signal term of the Wiener filter 
is marginalized over the power spectra constrained by the data of 
this patch themselves (Eriksen et al. 2004; Larson et al. 2007). The 
maps include only modes for I > 76 and smoothed to 1° FWHM; 
lower multipoles are removed due to a significant foreground con- 
tribution. Note the clear difference in amplitude: the E modes 
show a high-signal-to-noise cosmological signal while the B modes 
are consistent with noise. Maps for the other patches are available 
online. 



K-band (23-GHz) map through pipeline A and estimate 
its band power, C^ K , as well as the cross spectra with 

the QUIET Q-band data, d\ QK . These results are shown 
for the first bin (25 < I < 75; 6=1) in Table 6, 
together with the corresponding QUIET band powers, 
Cjf . Since foregrounds do not contribute to the sample 

variance, the uncertainties for and C*^^ are given 

by instrumental noise only, including contributions from 
both WMAP and QUIET. For C^, sample variance as 
predicted by the ACDM model is also included. 

There is significant EE power in patch CMB-1 as mea- 
sured by C*^^. We also find a correspondingly signifi- 
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Fig. 11. — The top panel shows EE results with 68% C.L. error bars; the bottom panel shows BB 95% C.L. upper limits. For comparison, 
we also plot results from previous experiments (Brown et al. 2009; Chiang et al. 2010; Larson et al. 2010) and the ACDM model (the value 
r = 0.2 is currently the best 95% C.L. limit on tensor modes). 



cant cross correlation between the WMAP K band and 
the QUIET Q band, confirming that this excess power is 
not due to systematic effects in either experiment and is 
very likely a foreground. No significant power is found in 
any other case. The non-detection of foreground power 
at I > 75 is consistent with the expected foreground de- 
pendence: cx £~ 2 - 5 (Carretti et al. 2010), and the low 
power found in C*^. 

The excess power observed in the first EE bin of CMB- 
1 is fully consistent with a typical synchrotron frequency 
spectrum. To see this, we extrapolate from K band 
to Q band, assuming a spectral index of (3 = —3.1 (Dunk- 
ley et al. 2009), and calculate the expected power in 



and C, 



QQ 

6=1' 



QK = U05 / 43.1 
b=1 1.01 V 23 



= 2.57 ± 0.69 (iX z 



(8) 



C, 



QQ 

6=1 



1.05 / 43.1 
L01 \~W 



Cf=i = 0.38 ± 0.10 fiK z , (9) 



where the prefactor accounts for the fact that (3 is de- 
fined in units of antenna temperature, and the uncer- 
tainties are scaled from that of C™. These predictions 
are fully consistent with the observed values of and 

C*6=d' wnen combined with the ACDM-expected power. 
We conclude that the excess power is indeed due to syn- 
chrotron emission. 

7.3. Constraints on Primordial B modes 

We constrain the tensor-to-scalar ratio, r, using the 
QUIET measurement of the BB power spectrum at low 
multipolcs (25 < £ < 175). Here r is defined as the 
ratio of the primordial-gravitational-wave amplitude to 
the curvature-perturbation amplitude at a scale ko = 
0.002 Mpc -1 . We then fit our measurement to a BB- 
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Fig. 12. — CMB power spectra are shown for each patch individ- 
ually. The top and bottom panels show the EE and BB spectra, 
respectively. The different error bars for each patch mainly reflect 
the amounts of time each was observed. 



For pipeline A, we find r = 0.35lg'gy, correspond- 
ing to r < 2.2 at 95% confidence. Pipeline B obtains 
r = 0.52^q;|^. The results are consistent; the lower panel 
of Figure 11 shows our limits on BB power in comparison 
with those from BICEP, QUaD, and WMAP. QUIET lies 
between BICEP and WMAP in significantly limiting r 
from measurements of CMB-B-mode power in our mul- 
tipole range. Although we neither expected nor detected 
any BB foreground power, the detection of an EE fore- 
ground in patch CMB-1 suggests that BB foregrounds 
might be present at a smaller level. We emphasize that 
the upper limit we report is therefore conservative. 

7.4. Temperature Power Spectra 

Figure 13 compares the QUIET and WMAP Q-band 
temperature maps and TT, TE, and TB power spectra. 
Agreement with the ACDM model is good. This is a 
strong demonstration of the raw sensitivity of the QUIET 
detectors; the single QUIET differential-temperature as- 
sembly produces a high-signal-to-noise map using only 
189 hours (after selection) of observations. The high sen- 
sitivity of these modules makes them very useful for cali- 
bration, pointing estimation, and consistency checks (see 
QUIET T WMAP T 
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TABLE 6 
Band and Cross Powers for 



25-75 



Patch 


Spectrum 


/.KK 
°6=1 




u b=i 


CMB-1 


EE 
1313 
EB 


17.4 ± 4.7 

4.8 ±4.5 
-6.2 ±3.2 


3.30 ±0.55 

0.40 ±0.41 
0.27 ±0.38 


0.55 ±0.14 

0.06 ± 0.08 
0.10 ± 0.08 


CMB-2 


EE 
BB 
EB 


5.5 ±3.7 

4.6 ±3.4 
-5.5 ±2.8 


0.01 ±0.56 
0.18 ±0.48 
-0.39 ±0.41 


0.23 ±0.19 
-0.11 ± 0.13 
-0.20 ±0.12 


CMB-3 


EE 
1313 
EB 


0.2 ±1.9 
-0.3 ±2.6 
1.4 ± 1.7 


0.64 ±0.43 
0.33 ±0.35 
-0.34 ±0.30 


0.10 ±0.18 
0.01 ± 0.13 
-0.27± 0.11 


CMB-4 


EE 
1313 
EB 


-5.2 ±5.1 
-2.6 ±5.2 
-1.0 ±3.9 


0.7 ± 1.2 
-0.1 ± 1.1 
0.0 ±0.9 


0.65 ± 0.58 
-0.37 ±0.52 
-0.15 ± 0.47 



Note. — Power-spectra estimates for the first multipole bin 
for each patch, computed from the WMAP7 K-band data and the 
QUIET Q-band data. The units are 1(1 + l)C f /27r (a*K 2 ) in ther- 
modynamic temperature. Uncertainties for C^J^ and C^J^ include 

noise only. For C^S^ they additionally include CMB sample vari- 
ance as predicted by ACDM. Values in bold are more than 2 a away 
from zero. 

spectrum template computed from the ACDM concor- 
dance parameters with r allowed to vary. For simplicity, 
we fix the tensor spectral index at nt — in comput- 
ing the template 42 . This choice makes the BB-power- 
spectrum amplitude directly proportional to r. 

42 Our definition of r agrees with Chiang ct al. (2010) 
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Fig. 13. — The top row compares our temperature map to the 
WMAP 7-year Q-band map (Jarosik et al. 2010) for patch CMB-1 
in Galactic coordinates. Lower panels show the CMB temperature 
power spectra: TT, TE, and TB. 
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Section 4). 



8. SYSTEMATIC ERRORS 

The passing of the null suite itself limits systematic 
uncertainty, but to get well below the statistical errors, 
dedicated studies are needed. They are important in 
gaining confidence in the result and also in evaluating 
the potential of the methods and techniques we use for 
future efforts. We pay special attention to effects that 
can generate false B-mode signals. Our methodology is 
to simulate and then propagate calibration uncertainties 
(see Section 4) and other systematic effects through the 
entire pipeline. The systematic errors in the power spec- 
tra are shown in Figure 14. The possible contaminations 
are well below the statistical errors; in particular, the 
levels of spurious B modes are less than the signal of 
r = 0.1. This is the lowest level of BB contamination yet 
reported by any CMB experiment. This section describes 
how each effect in Figure 14 is determined and considers 
three additional possible sources of contamination. 

An uncertainty not shown in Figure 14 is that aris- 
ing from the overall responsivity error estimate of 6% 
(12% in power-spectra units). After including the effect 
of possible time-dependent responsivity variations (4%, 
see below), the power-spectra uncertainty is 13%. It is 
multiplicative, affecting all power-spectra results inde- 
pendent of multipole. 

8.1. Beam Window Function and Pointing 

The uncertainty in the beam window function is an- 
other multiplicative factor, one which increases with mul- 
tipole. We estimate this uncertainty using the difference 
of the beam window functions measured for the central 
module and the modules of the differential-temperature 
assembly, which are at the edge of the array. The differ- 
ence is statistically significant, coming from the different 
locations (with respect to the optics) in the focal plane; 
it is expected from the pre-season antenna range mea- 
surements. 

Uncertainties in pointing lead to distortions in polar- 
ization maps. E power will be underestimated and spuri- 
ous B power (if the distortions are non-linear) generated 
(Hu et al. 2003). We quantify these effects by using the 
differences in pointing solutions from two independent 
models: the fiducial model used for the analysis and an 
alternative model based on a different set of calibrating 
observations. We also modeled and included the effects 
of the deck-angle-encoder shift which occurred for a por- 
tion of the season (Section 4.3). 



8.2. Responsivity and Polarization Angle 

Responsivity shifts, particularly within CESes, lead to 
distortions in the maps. Full-pipeline simulations quan- 
tify the shifts caused by variations in the cryostat or 
electronics temperatures. Similarly shifts from using re- 
sponsivities determined from the Moon data, Tau A data, 
or from the sparse wire grid, rather than those from the 
sky dips, are determined. We also incorporate the un- 
certainty in the atmospheric-temperature model used in 



analyzing the sky-dip data. The largest possible effects 
on the power spectra are shown in Figure 14. 

Uncertainties in the orientation of the polarization axes 
of the modules can lead to leakage between E and B 
modes. To quantify this leakage, we use the differences 
in power spectra where these angles are determined from 
Moon data, Tau A data, and the sparse-wire-grid data. 
As expected, the largest effects show up in EB power. We 
also estimate systematic error due to possible fluctuation 
of the detector angles over the course of the season. The 
contribution is negligibly small compared to the overall 
shift of the angles described above. Both effects are in- 
cluded in the "Polarization Angle" points in Figure 14. 

8.3. Instrumental Polarization 

As described in Section 4.5, the I to Q (U) leakage 
coefficients for the QUIET detector diodes are small: 
1% (0.2%). Except in the case of patch CMB-4, our 
scanning strategy significantly reduces this effect with 
the combination of sky and deck-angle rotation. 

We estimate spurious Q and U in the maps for each 
CES-diode using the WMAP temperature map and our 
known leakages. Shown in Figure 14 are the estimates of 
spurious EE, BB, and EB powers from full-pipeline sim- 
ulations, where for each realization the spurious Q and 
U are added to the Q and U from simulated ACDM E 
modes. While this method has an advantage of being 
able to use the real (not simulated) temperature map, it 
does not incorporate TE correlation, which only affects 
the spurious EE power. As a complement, we repeat the 
study, but using simulated ACDM maps for both tem- 
perature and polarization; this only changes the estimate 
of spurious EE power by 30% at most. Because the spu- 
rious power is as small as it is, we have treated it as a 
systematic rather than correcting for it. Doing so would 
give us a further order of magnitude suppression. 

Differing beam ellipticities can also induce higher mul- 
tipole polarization signals. We measure these leakages 
from Tau A and Jupiter observations and find that the 
higher-order multipolcs arc at most 0.1% of the main- 
beam peak amplitude. The corresponding effects on the 
power spectra, which are seen in Figure 14, are of little 
concern. 

8.4. Far Sidelobes Seeing the Sun 

While we make cuts to reduce the effects of far side- 
lobes seeing the Sun (Sections 2 and 5.1.3), small con- 
taminations could remain. We make full-season maps 
for each diode in Sun-centered coordinates and then use 
these maps to add contamination to full-pipeline CMB 
simulations. The excess power found in the simulations 
is taken as the systematic uncertainty. We do not ob- 
serve any signature from the Moon, nor do we cut on 
proximity to it. We estimate the related systematic er- 
ror and find that it is negligibly small compared to that 
assigned to the contamination from the Sun. 

8.5. Other Possible Sources of Systematic Uncertainty 

Here we discuss a few additional potential sources of 
systematic uncertainty, which are found to be subdomi- 
nant. 

Ground-Synchronous Signals. QUIET's far side- 
lobes do see the ground for some diodes at particular 
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Fig. 14. — Systematic uncertainty estimates for EE, BB, and EB power spectra. Estimates for a variety of effects (see text) are shown 
for the three power spectra. In all cases, they are well below the statistical errors, which are also shown. In particular, the contaminations 
to the primordial— B-mode signal, at multipoles below 100, are below the level of r = 0.1, even though we do not make a correction for the 
largest contaminant, the monopole leakage. 



elevations and deck angles. Ground pickup that is con- 
stant throughout a CES is removed by our TOD filters; 
the net effect of this filtering in the full-season maps is a 
correction of sa 1 fiK. 

The only concern is ground pickup that changes over 
the short span of a single CES. We find little evidence for 
changes even over the entire season, let alone over a single 
CES. We therefore conservatively place an upper limit on 
such changes using the statistical errors on the ground- 
synchronous signal. We start with the CES and module 
with the largest ground pickup. We then simulate one 
day's worth of data, inserting a ground-synchronous sig- 
nal that changes by its statistical error. Given the dis- 
tribution in the magnitude of the ground-synchronous 
signal and assuming that changes in this signal are pro- 
portional to the size of the signal itself, by considering 
that the signals from changing pickup add incoherently 
into the maps made from multiple CES-diodcs at a vari- 
ety of elevations and deck angles, we estimate an upper 
limit on residual B power from possible changing ground- 
pickup signals. The result is < 10~ 4 /iK 2 at multipoles 
below 100. 

ADC Non-linearities. The possible residual after 
the correction for the non-linearity in the ADC system 
results in effects similar to the I to Q (or U) leakage 
and the variation of the responsivity during the CES. 
We estimate such effects based on the uncertainty in the 
correction parameters, confirming that there is at most 
a 3% additional effect for the leakage bias, and that the 
responsivity effect is also small, less than half of the sys- 
tematic error shown for the responsivity in Figure 14. 

Data-Selection Biases. Cuts can cause biases if they 
are, for example, too stringent. We expect none but to 
be sure we apply our selection criteria to 144 CMB + 



noise simulations. No bias is seen, and in particular we 
limit any possible spurious B modes from this source to 
< 10 -3 yK 2 at multipoles below 100. 

9. CONCLUSIONS 

QUIET detects polarization in the EE power spectrum 
at 43 GHz. We confirm with high significance the detec- 
tion of polarization in the region of the first acoustic peak 
(Chiang et al. 2010) in the multipole region i = 76-175. 
We find no significant power in either BB or EB between 
£ = 25 and I = 475. We measure the tensor-to-scalar 
ratio to be r = 0.35±J^7- 

These results are supported by a very extensive suite of 
null tests in which 42 divisions of data were used for each 
of 33 different cut configurations. The selection criteria 
and systematic errors were determined before the power 
spectra themselves were examined. Biases were revealed 
during this process, the last of which was a contami- 
nation present in the null spectra at the level of about 
20% of the statistical errors, but eliminated when cross- 
correlating maps with differing telescope pointings. The 
robustness of the final results is further supported by 
having two pipelines with results in excellent agreement, 
even though one uses only cross correlations while the 
other also uses auto correlations. 

Several possible systematic effects are studied with full 
end-to-end simulations. The possible contaminations in 
the B-mode power are thereby limited to a level smaller 
than for any other published experiment: below the level 
of r = 0.1 for the primordial B modes; simply cor- 
recting for the known level of instrumental polarization 
would reduce this to r < 0.03. This very low level of 
systematic uncertainty comes from the combination of 
several important design features, including a new time- 



18 



The QUIET Collaboration 



stream "double-demodulation" technique, side-fed Drag- 
onian optics, natural sky rotation, and frequent deck ro- 
tation. 

The correlation modules we use have a polarization 
sensitivity (Q and U combined) of 280 /iK-y/s, leading to 
an array sensitivity of 69 ^K-y/s. Further, the 1// noise 
observed in our detectors is small: the median knee fre- 
quency is just 5.5 mHz. One important outcome of this 
work, then, is the demonstration that our detectors, ob- 
serving from a mid-latitude site, give excellent sensitivity 
and systematic immunity. 

Because of our mid-latitude site, we are driven to col- 
lect data in four separate patches. While we lose some 
sensitivity (compared to going deeper on a single patch), 
there are a few advantages that we have exploited. The 
patches are scanned differently, in terms of time of day 
and the degree of crosslinking, and these differences allow 
some important systematic checks. Another advantage 
concerns foregrounds. 

Foreground contamination is expected to be one of 
the main limiting factors in the search for primordial B 
modes. Indeed we report a 3-ct detection of synchrotron 
emission in one of our four CMB patches, originally cho- 
sen for their expected low foreground levels. Our detec- 
tion is only in EE; our BB 2-er limit permits a BB signal 
about half as large. If we extrapolate that value to the 
foreground minimum of about 95 GHz, we would have 
synchrotron contamination at the level of r = 0.05 from 
this one patch. Neither WMAP nor Planck will have 
enough sensitivity (Tauber et al. 2010) to sufficiently con- 
strain the polarized synchrotron amplitude at this level 
in any patch. In fact, our Q-band polarization maps arc 
already as deep or deeper than what Planck will achieve 
at the same frequency. Dedicated low-frequency obser- 
vations are clearly needed to achieve such constraints. 
When foreground cleaning becomes important, consis- 
tency among separate patches will be an important han- 
dle on our understanding. 

Further progress must be made through larger arrays 
and longer integration times. In hand we have data col- 
lected by the 90-element W-band array with similar sen- 
sitivity to our Q-band array and more than twice the 
number of observing hours. Results from the analysis of 
that data set will be reported in future publications. A 
W-band receiver with the sensitivity to reach below the 
level of r = 0.01 is under development. 
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the Strategic Alliance for the Implementation of New 
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